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Abstract 

We investigate the kinetics of systems in which particles of one species undergo 
binary fragmentation and pair annihilation. In the latter, non-linear process, frag- 
ments react at collision to produce an inert species, causing loss of mass. We analyse 
these systems in the reaction-limited regime by solving a continuous model within 
the mean-field approximation. The rate of fragmentation, for a particle of mass x 
to break into fragments of masses y and x — y, has the form x A_1 (A > 0), and the 
annihilation rate is constant and independent of the masses of the reactants. We 
find that the asymptotic regime is characterized by the annihilation of small-mass 
clusters, with the cluster-density decaying as in pure-annihilation and the average 
cluster-mass as in pure-fragmentation. The results are compared with those for a 
model with linear mass-loss (i.e. with a sink rather than a reaction). 

We also study more complex models, in which the processes of fragmentation and 
annihilation are controlled by mutually-reacting catalysts. Both pair- and linear- 
annihilation are considered. Depending on the specific model and initial densities of 
the catalysts, the time-decay of the cluster-density can now be very unconventional 
and even non-universal. The interplay between the fragmentation and annihilation 
processes and the existence of a scaling regime are determined by the asymptotic 
behaviour of the average-mass and of the mass-density, which may either decay 
indefinitely or tend to a constant value. We discuss further developments of this 
class of models and their potential applications. 
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1 Introduction 



Irreversible reaction-diffusion processes, such as aggregation or annihilation of species, 
occur in many physical, chemical and biological systems. In the past decade they 
have been investigated extensively. This followed from the realization that spatial 
fluctuations in the reactant densities slow down the kinetics and lead to an asymp- 
totic behaviour at low- dimensions which disagrees with the classical reaction-rate 
equations. (The upper critical dimension, d c , is defined as the spatial dimension be- 
low which such disagreement occurs.) Recently, a combined process of aggregation- 
annihilation of two species has been studied exactly in 1 dimension by Krapivsky 
DH, and numerically in higher dimensions by Sokolov and Blumen |f|. A class of 
related models, but with an arbitrary number of species, has also been investigated 
at the mean-field level by Ben-Nairn and Krapivsky Jl| . The power-law decay of the 
cluster-densities has been understood, and in many cases the cluster-mass distribu- 
tions have been found to approach universal scaling forms. 

Fragmentation is another irreversible kinetic process with important applications 
in nature and technology. These include: polymer degradation, through shear action 
HI, chemical attack ||, or exposure to radiation ||; droplet breakup |7j; fibre- 
length reduction ||; fragmentation of colloidal aggregates and rock crushing 
and grinding M. New insight into the statistical mechanics of fragmentation has 



followed from the work of Ziff and McGrady [[U], [II], who solved exactly a class of 
models with different fragmentation kernels. 

The aggregation-fragmentation (or reversible-aggregation) process has been con- 
sidered by many authors. This has been applied to studying reversible polymer- 
ization, for instance, by van Dongen and Ernst [l2j and by Sintes, Toral and 
Chakrabarti |jl3| . More generally, mean-field approaches have been used to study 
the steady-state scaling properties of reversible-aggregation, for example, by Family, 
Meakin and Deutch |14| and by Blackman and Marshall In both aggregation- 
annihilation and aggregation-fragmentation processes, rich kinetic behaviours arise 
due to the competing reaction mechanisms. 

In this paper we focus on fragmentation processes with mass-loss. A number 
of authors have considered models combining fragmentation and linear- annihilation 
processes (i.e. a sink, where a species 'evaporates' at a constant rate in time) (see, 
e.g., [JTSL [T7|, |18[ and references therein). These models are relevant, for example, 
to studying oxidation and dissolution of solid porous media |I7|, [OJ. However, to 
the best of our knowledge, the interplay between annihilation and fragmentation re- 
mains unexplored in the case of binary- annihilation reactions. It is this subject with 
which we are concerned. The kinetics of these phenomena are non-linear, and there- 
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fore more complex and, in principle, more difficult to analyse. Fragmenting systems 
where fragments react at collision (maybe in the presence of some agent) to produce 
an inert species, would yield potential applications of these combined processes. We 
investigate, at the mean-field level, a class of (non-linear) annihilation-fragmentation 
models, aimed at elucidating the competition between these processes. For all mod- 
els considered we include, for comparison, an analysis of the corresponding model 
with linear-annihilation (a sink). 

In the simplest situation of a single species, the cluster-density decays as in the 
case of pure-annihilation, while the cluster-mass distribution approaches an expo- 
nential scaling distribution. An exponential decay of the scaling distribution has also 
been found in pure-fragmentation, aggregation- annihilation |], § an d reversible- 
aggregation [14|. To allow for more interesting, and possibly more realistic be- 
haviour, we introduce additional species which work as catalysts for annihilation 
and/or fragmentation; reactions between the catalysts are also allowed for. We find 
that the density-decay of the main-species can be surprisingly unconventional, de- 
pending on the specific model and on the initial densities of the catalysts. The decay 
may be universal but contain logarithmic factors, or it may contain a power-law in 
time with a non-universal exponent which depends on the reaction rates. The inter- 
play between annihilation and fragmentation and the existence of a scaling regime, 
are determined by the asymptotic behaviour of the mass-density and typical cluster- 
mass. These may either decay indefinitely or tend to a constant value, depending 
on the relative amounts of annihilation and fragmentation catalysts. 

This paper is organized as follows. In Sec. 2 we introduce and solve a model 
with single-species fragmentation-annihilation where the fragmentation kernel has a 
rather general form. The results for the corresponding model with linear annihilation 
are stated for comparison. In Sec. 3 we include catalysts in our model and solve it for 
general time- dependences of the catalysts. Specific asymptotic results are presented 
in the case where the fragmentation and annihilation catalysts mutually annihilate. 
The conclusions and further comments and suggestions are given in Sec. 4. 



2 Single-species annihilation- fragmentation 
2.1 non-linear annihilation (model 1) 

The more elementary model combining the processes of pair-annihilation and frag- 
mentation is that involving a single 'molecular' species. This is described by the 
binary reaction scheme 
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where A4 denotes a cluster consisting of i monomers of species A. The following 
calculations could be repeated for a discrete model, but for analytical convenience 



we shall take the continuous limit, which is standard practice pA\ , [19]. 

Let a(x, t) be the density distribution of cluster-masses (or sizes) x at time t, 
and 

POO 

A n (t) = / dxx n a(x,t) (2) 
Jo 

its moments. In particular, the zeroth and first moments, A(t) = Ao(t) and Ai(t), 
denote the total number- density of clusters and the total mass-density, respectively. 
For simplicity, we consider linear fragmentation (driven by a homogeneous external 
agent [|T^]), neglect multiple breakup, and take the rate that a cluster of mass x 
breaks into clusters of mass x — y and y, F(y,x — y), to be a constant (which we 
set to be unity by a rescaling of time). Hence, the total rate of fragmentation of 
any cluster, Jq dyF(y, x — y), is equal to its mass, x. We also envisage the simplest 
possible reaction, where two clusters of any mass mutually annihilate upon encounter 
at a rate independent of their masses. Later, we shall see how to generalize our 
results for a kernel of the form F(y, x — y) = x x ~ l , independent of the masses of the 
fragments. 

Treating the annihilation term within the mean-field description, which is valid 
in the reaction-controlled regime, the kinetic equation for a(x, t) reads 

POO 

d(x) = —xa(x) + 2 / dya(y) — Ja(x)A , (3) 



where J/2 is the relative rate of annihilation (equal rates of fragmentation and 
annihilation correspond to J = 2). A set of non-linear equations for the moments 
of a(x) follows immediately, 

1 — n 

A n = -—A n+1 - JA n A . (4) 
1 + n 

A closed set of equations is found for the first two moments, 

A = A X -JA 2 , (5) 
ii = -JAA X , (6) 

whilst the higher moments obey an infinite hierarchy of equations which are, ap- 
parently, not soluble. The absence of a linear term in ([|) is a reflection of the 
conservation of mass in the fragmentation process. Adopting more general rates 
of fragmentation and annihilation, say x x and (xy) a , respectively, would yield the 
following generalization of (^): 

A n = - — ; ~~ J A n + a A a . (7) 

1 + n 
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Clearly, no closed set of equations can be obtained for A and Ai, unless A = 1 and 
a = 0, which justifies our present choice. 

To solve for A(t) and Ai(t), we proceed as follows. From the cluster 

average-mass, L(t) = A\(t)/A(t), obeys the equation 



with solution 



L(t) 



6 + t 



(t 



oo 



(8) 



(9) 



where 5 = A /A w = 1/L(0) and A = A(0) and A w = A^O). Using this result, 
Eqs. @-@ can be solved to give 

A + A 1Q t 2 



A(t) 

Mt) 



1 + J[A t + A w t 2 /2] 
A{t) ^ 2 
J+t ~ ' 



Jt 



(10) 
(11) 



The asymptotic behaviour in ([L0| )-(11), and hereafter, hold at large t. 

To determine the mass distribution a(x,t), we map the problem onto a pure- 
fragmentation one using the auxiliary function 

rt 



b(x, t) = a(x, t) exp 



J [ dt'Ait') 
Jo 



The exponential factor is known and from (|10|) we see immediately that 



exp 



J /' dt'Ait') 
Jo 



J[A t + A 10 t 2 /2] . 



(12) 



(13) 



The function b(x, t) obeys a standard fragmentation equation, i.e. it is the same as 
@ but without the non-linear annihilation term. Its general solution is well known 
0, p0[, and can be written as 



b(x,t) 



b(x,0) + y°° dyb(y,0)[2t + t 2 (y - x)]^j exp(-te) . 



(14) 



We find it easier and more illustrative, however, to work from a particular initial 
condition, and then to extract the asymptotic, universal behaviour from the special 
solution obtained. We adopt the exponential initial distribution 



a(x, 0) = b(x, 0) = Pq exp(— box) 



(15) 



where b = S and (3q = A 6. The corresponding solution is obtained by substituting 
(|15D in (jl4l) or, equivalently, attempting the ansatz b(x,t) = (3{t) exp[— b(t)x] in the 
equation for b(x,t) |pl| . This gives, after including the time-dependent factor (^3|), 

A w (5 + t) 2 . (s , . . A(t) 



a(x, t) 



J[A t + A 10 P/2] 



exp [—(5 + t)x] 



Lit) 



exp[-(5 + t)x] . (16) 
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As t — > oo, a(x,t) ~ (2/ J) exp(— tx) becomes independent of the initial condition. 
This universal asymptotic behaviour can also be obtained by looking for a scaling 
solution of the form 

a(x,t) = L-^(C) , ( = x/L, (17) 

which we expect to hold for general initial conditions, in the regime when t is large 
and x is small, with xt an arbitrary constant. From (JlQ|)-(Jll|) at large t, we find 
w — 0. Replacing a by in Eq. (|3|) and taking the limit t — > oo, yields a linear 
differential equation with solution 

a(x,t) = 0(C) = jexp(-C) • (18) 

One can calculate all the moments of a(x,t) given by the special solution (|16f). 
This yields for and Al(£) exactly the same expressions as in (0)-(0), which 
hold for general initial conditions (One can check that this is also true for modified 
exponential initial conditions, e.g. x p exp[— b x].). Since Eq. (|4]) for all moments is 
recursive, we conclude that the special solution ([16]) must have the same moments 
as the general solution. Hence, we did not need to solve for A(t) in the first place; 
rather, we could have solved self-consistently for A(t) from this special distribution. 
The rationale is: that one may be able to determine a particular mass distribution 
even when A(t) is not known a priori, i.e. even when there is not a closed set of 
equations as in the more general case of Eq. (|7|). 

With this insight, we shall now solve the problem for a more general fragmenta- 
tion kernel, F(y,x — y) — x A-1 , with A > 0. There is no reason why the previous 
statement should not hold here, i.e. we expect that for general A a particular expo- 
nential solution will have the same moments as the general solution. Eq. (|3]) now 
becomes 

a{x) = -x x a(x) + 2 dyy x ~ x a{y) - Ja(x)A . (19) 

J X 

Once again we use the auxiliary function b(x,t), given by (|12"D, although now the 
exponential factor is not known. This gives a pure- fragmentation equation for b(x, t), 

b( x ) = -x x b{x) + 2 / dyy x - x b{y) , (20) 

J X 

which has a known, but complicated, general solution [l(|. It has been proven that 
A > is a necessary and sufficient condition for b(x, t) to have a scaling, universal 
behaviour fl9], ^2f . We are allowed, therefore, to deduce this behaviour from a special 
solution to ([19]). A simple solution with exponential form (whose initial condition 
is an obvious generalization of ( fjjD ), is 

Kx, t) = ( 6o + t f' x exp[-(6 + t)x x ] , (21) 
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with moments 

Bn(t) = /^Vfm • (22) 
w (&o + t) (n ~ 1)/A 

In (pip and ( p2j ) appropriate normalization has been used. We now have 

&S /A = 5A 2) i , /3o = , (23) 

with the convenient notation 

A M = A m (X) ee . (24) 

To find the explicit solution to the problem, we must determine the density of 
clusters, A(t), self-consistently. From (p2|) and the transformation (|l|) we obtain 
an integral equation. This can be recast as the following differential equation 



with solution 



A(l + t/6 ) 1/A 1 + A 1 



m = i + jA^&Ki+t/W+w-i] * —Tt ■ ^°°> ■ (26) 

It is easy to see that exp(J J^dfA) is simply the denominator of A in fl26|) . From 
this, (O) and (E2), we find 



A + 1 1 

= A B+1 , 1(6o + t)B/A ~ - Jtl+n/A , (27) 

, . AA 2 i A(t) . , Al exp[-C A ] 

a 0M) = F7TTTTT7IV exp[-(6 + t)x Al J 



T(l/A) L{t) v u 7 J JtL(t) ' 

with C = h.2,1 x/L(t). All results reduce to the previous ones if we set A = 1. Apart 
from numerical factors, the essential modifications in the asymptotics can be alluded 
to the time-dependence of the average cluster-size, L ~ t _1 / A , and the dependence 
of a(x, t) on the scaling variable (. We see that an increase in A slows down the 
decay of L. This follows from the resulting decrease in the rate of fragmentation of 
small fragments. 

The mass distribution in fl2"T|), or (|T8|), is peaked around ( = 0, showing that 
the majority of the clusters have vanishingly- small mass at late times due to frag- 
mentation. The factor 1/J indicates that if the annihilation rate increases the 
mass and cluster densities should decrease. The cluster density A(t) is the only 
moment whose decay is independent of A. We recall that in pure-annihilation pro- 
cesses A ~ 1/Jt, while the cluster-mass distribution remains invariant in time, 
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apart from an overall time-dependent factor. Our results indicate, therefore, that 
the fragmentation-annihilation process at late-times consists essentially of the an- 
nihilation of small-mass clusters. This leads us to conjecture that in this case the 
upper critical dimension is 2, as in the pure-annihilating case. That is clearly true 
for a discrete system, where there is a minimum cluster-mass and the fragmentation 
process ends before all particles annihilate. 

2.2 linear annihilation (model 2) 

Here, we briefly discuss model 2: a model with the same fragmentation process as 
model 1 but with linear annihilation. In this case there are no binary reactions 
and single clusters annihilate at a constant rate in time. We introduce this model 
for comparison with model 1 and because in Sec. 3 we shall make an interesting 
generalization of it. 

The kinetic equation for the cluster- mass distribution is similar to (|19|), but 
with the annihilation term replaced by the linear term —Ja(x). This yields a trivial 
modification of the pure- fragmentation model ( pPf) for general A: the transformation 



([12]) is replaced by b(x,t) = a(x,t) exp[Jt], leading to a simple extra exponential 
factor in the results. Working once again with the special solution (pT|) , which 
corresponds to the initial condition a(x,0) = b(x, 0) = j3o exp[— box x ], we obtain 

A(t) = exp[-Jt] ~ t 1/x exp[-Jt] , (t -»• oo) , 

L{t) 

A n (t) = A n+lil ^-Al_ „ ^P_^ j (28) 

n(rrf\ AA 2,1 A (t) r , ,x Xl exp[~C A ] expJ-Jt] 

where L(t) is the same as in model 1 (Eq. (p7|)). Similar remarks about the depen- 
dence on A apply here. 



3 Annihilation-fragmentation with catalysts 

Model 1 has elucidated some of the features in the interplay between the fragmenta- 
tion and annihilation processes. However, it appears to have a somewhat unrealistic 
drawback: when two clusters of species A meet they react to produce a new, in- 
ert species, rather than merging into a larger cluster. Such a reaction would seem 
plausible, though, if another species, say C, was involved. In a reaction of the form 
A x + C — > inert, or A x + A y + C — > inert, the kinetic equations of A and C would 
be coupled together. To be able to make progress, however, we envisage a situation 
where the reactions are induced by the presence of C but C does not take part in 
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the reactions, i.e. C is a catalyst. In Sec. 2, it was implicitly assumed that annihi- 
lation was driven by a catalyst of constant concentration. Here we shall allow the 
concentration of C to vary in time. Catalysts play a key role in chemical, biological 
and ecological systems (see, e.g., [^, [24j]), so the study of this class of models is 
potentially of great interest. 

With these insights, let us consider the following two models where the fragmen- 
tation of A is controlled by a catalyst B, 

A x+y + B^A X + A y + B, (29) 

and the annihilation of A is induced by a catalyst C. With a non-linear annihilation 
reaction (model CI) 

A x + A y + C -> inert + C , (30) 
while with a linear annihilation reaction (model C2) 

A x + C -> inert + C . (31) 

In addition, we allow the catalysts to undergo the following annihilation reaction, 
at rate R: 

B + inert . (32) 

In the reaction-controlled regime B(t) and C(t), the concentrations of species B and 
C, have the following time dependence (see, e.g., PSJ): 



B(t) = C(t) = l + B ° aRt , (B„ = C„) , 

B(t) = AB ^P(- A «) , C(t) = A + B(t) , (Co>Bo) , (33) 
C — B exp(-ARt) 

with B = B(0), C = C(0) and A = \C - B \. Note that if B(t) and C(t) were 
constant in time, then, via a trivial rescaling, models CI and C2 would reduce to 
models 1 and 2 without catalysts. 

3.1 method of solution 

We shall solve the generalizations of models 1 and 2, with general values of A, 
following a similar procedure to that used in Sec. 2. The essential difference now, is 
that a formal time variable has to be defined, r(t), which plays the role of time t in 
the previous models. 

Model CI, for non-linear annihilation in the presence of catalysts, is defined 
by the reaction schemes (^) and (p0[). The kinetic equation for the cluster-mass 
distribution, which generalizes flT9]), is 



a(x) 



x x a(x) + 2 / dyy x ~ l a(y) 



B - Ja(x)AC , (34) 



8 



while the equation for the moments of a(x, t) now reads (cf. 

1 — n 



with a = 0) 



A, 



1 + n 



Ai+a B - JA n AC 



(35) 



For A = 1, a closed set of equations can be obtained once again, from which the 
general solution for A(t) and Ai(t) can be found. For general A, we map ([54]) onto 
a pure-fragmentation equation using the transformation (cf. (|T2"|)) 



b(x, t) = a(x, t) exp 
and redefining time as 

T(t) = 



J I dt'A(t')C(t') 
Jo 



(36) 



dt'B(t') 



(37) 

Clearly, one only expects a scaling regime to occur if r(£) — > oo as £ — > oo. Then 
r) obeys Eq. (p0|). As usual, we work with the special exponential solution ( pip 



(now with r replacing £). Inserting this solution into (36) and differentiating once, 
yields the self-consistent equation (cf. (p5f)-(p6|)) 

AB 



with solution 



A(t) 



A — — — r J A C , 

A(6 + r) 
AoAi^o + r) 1 ^ 



1 + JA 10 A h2 f* dt'C(t') % + r(f )] V> • 
As in Sec. 2, exp(J J*q dt'AC) is simply the denominator of A in 
(H), (H) and © we find 

A 2 ,i 



(38) 
(39) 

Hence, from 



Lit) 

A n {t) 

a(x, t) 



A(t) 



A 



"+ 1 ' 1 ( 6o + T )n/A 

AA 21 A{t) 



(40) 



exp[-(6 + r)a; ; 



T(l/A) Lit) 

The essential differences from ( PTD are due to the presence of r(t). 

We now turn to model C2, for linear annihilation in the presence of catalysts, 
defined by the reaction schemes ( |29| ) and (|3T|). The kinetic equation for a(x,t) is 

similar to (|54]), but with the non-linear term replaced by — Ja(x)C. Defining 

ft 



J [ dt'Cit') 
Jo 



(41) 



b(x, t) = a(x, t) exp 

and using the same formal time variable r(t), b(x,r) obeys Eq. (|20|) once again, 
and we work with the special exponential solution as before. The results, therefore, 
follow from (pl|)-(p2|) (with r replacing t) multiplied by exp[— J f dt'C]. This gives 

(cf. (H) 

Ait) = 4^ exp[-J ['dt'Cit')} , (42) 
' ' ' Jo 



Lit) 



with L(t), A n (t) and a(x,t) still given by (|4Q|) . 
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3.2 asymptotic results for B + C — > meri 

The results derived in Sec. 3.1 for models CI and C2 hold for any time dependence 
of the catalysts. To obtain the asymptotic behaviours corresponding to catalysts B 
and C mutually annihilating according to (^), we insert the mean-field expressions 
Q3"3|) into the previous results and look for the dominant contributions in the limit 
t — > oo. There are three qualitatively different cases, depending on the relative 
values of the initial densities of the catalysts. 



Bp = C 

With balanced amounts of the two catalysts, we have 



r = ln(l + RB t) 1/R ~ ln(t 1/R ) , (t -> oo) . (43) 



This gives, 



A 2 ,ifl 1/A 
[ln(t)] 1 

A n (t) ~ L(0"^(0 , (44) 
Aft) 

a M) ~ T7r(exp[-C A ] , C = V/L(f), 



L(t) 



for models CI and C2, and 



~ I±^_^ ^ 1 , model CI 



A Jln(t) ln{t J / R ) 
Aig 1 pn(t)]V* 
L(t) [5 #£] J/i? ~ t J /« 



model C2 . (45) 



As in models 1 and 2 without catalysts, the late-time regime is dominated by the 
annihilation of small-mass clusters and there is scaling. The cluster-number den- 
sity, however, decays logarithmically in model CI and with a non-universal power- 
law in model C2: the non-universal exponent depends on the annihilation rates 
J and R of species A and of the catalysts. We conclude that, although different 
time-dependences emerge, the balance between the annihilation and fragmentation 
processes in this symmetric case is analogous to that occurring in the absence of 
catalysts. 

Bp > C 

If there is a majority of the fragmentation catalyst the auxiliary time r behaves 
essentially like t, 

t = At + \n{(B -C exp[-ARt])/A} 1/R ~ At , (t -»• oo) . (46) 
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The key point to notice when evaluating the integral in the denominator of (^) and 
the exponential term in (42), is that J c?i / C(t / ) — ► ln(I?o/A) 1//i? as t — > oo. This 



yields, for models CI and C2, 



where 



L(t) " (^A ' (* ^ °°) ' 

A(t) ~ A 10 F(A)A 2il (At)V> , (47) 

- A 10 F(A) , 

~ jjj^ exp[-C A ] , ( = A 2)1 x/L(t) , 



F ^'> = i + jaU./- • m ° del 01 

A v J/R 

model C2 (4* 



1 oo 



Bo 

/ b + x + \n{C /(B - Aexp\Rx})} 1/R . (49) 

jo L J 

In this case the fragmentation process is much faster than the annihilation one, 
which is suppressed as the mass-density tends to a fraction F(X) of its initial value 
Aiq. Once again there a scaling regime. In conclusion, the asymptotic behaviours of 
models CI and C2 are similar to that of a pure-fragmentation process with rescaled 
initial mass-density A 10 F(X) and rescaled time At, where A is the residue of cata- 
lyst B. Hence, the smaller the value of A the longer the time needed to reach this 
regime. The value of F(X) is determined not only by the kinetic parameters J, R 
and A, but also by the initial densities B Q and C . 

Cp > B 

With a majority of the annihilation catalyst r tends to a finite value, 
r = lnUCo-^oexphAi^D/A} 1 ^ -> = ln(3,/A) 1/fl , (t -+ oo) . (50) 

Noticing that the integral in the denominator of (|39|), f dt'C(t')[b + r(t')] 1 / A , is 
dominated the large-t contribution, we obtain 

L(t) — ^ = A 2 ,i/(6o + ^oo) 1/A , (t -»• oo) , 
A n (t) ~ A(t) , (51) 
a(x,t) ~ A(t) exp[-(6 + Too)x A ] , 



for models CI and C2, and (cf. © and @8| 

1 

J At ' 

A 10 /A\ J / fi 



A(t) ~ , model CI 

v ; J At 



Loo \Cq 



exp[-JAt] , model C2 . (52) 
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Here, the annihilation process dominates over fragmentation. The latter is sup- 
pressed as the average cluster-mass tends to a non-zero value and the mass 
distribution becomes stationary (up to an overall time-dependent factor) and iden- 
tical to its initial form a(x, 0). This form of a(x, t) follows from all fragments being 
equally-likely to annihilate, and the prefactor A(t) — > indicates the eventual loss 
of all mass. This situation differs from B Q = Co, where annihilation occurs between 
increasingly smaller clusters. As expected, there is no scaling regime since r has a 
finite limit. In addition, the behaviour of a(x,t) may not be universal in this case. 
In conclusion, the asymptotic behaviour of models CI and C2 is similar to that of 
a pure-annihilation process with L = L m and rescaled time At: the process being a 
binary reaction for model CI, and linear with initial mass A w (A/Cq) j ^ r for model 
C2. The asymptotic average- mass, Loo, is determined by the kinetic parameters R 
and A, and by the initial values A /A w , B and Co- 

It should be stressed that the validity of the continuous-mass description used 
throughout is subject to restrictions in the long-time regime. In the cases Bq = Co 
and Cq > Bq, where the cluster- number A(t) — > and the system eventually reaches 
an empty state, the description holds while the system still contains a sufficiently 
large number of particles. This condition ensures the smoothness of the cluster-mass 
distribution. In the cases B = Co and B > Co, where the average cluster-mass 
L(t) — > and increasingly-smaller fragments are produced, the usual restrictions 
to a pure-fragmentation process apply [[II] . The essential difference within a dis- 



crete description lays in the behaviour of the density of the smaller fragments, or 
monomers. 

4 Conclusions 

We have solved a class of fragmentation models with non-linear mass-loss, at the 
mean-field level, and elucidated the competition between fragmentation and annihi- 
lation processes. Simultaneously, we have analysed the corresponding models with 
linear mass-loss. Although the solutions for the cluster-mass distribution hold for 
exponential initial conditions, we expect their moments to be identical to those of 
the corresponding general solutions. 

Models 1 and 2 (without catalysts) are characterized by the annihilation of small- 
mass clusters at late times. We believe that in this case the upper critical dimension 
is 2 (but see Q). 

Models CI and C2, with fragmentation and annihilation catalysts, have been 
solved for a general time-dependence of the catalyst densities. The linear annihila- 
tion model, in particular, shows much more interesting and less trivial behaviours in 
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this case. When B and C undergo mutual annihilation, we have found that differ- 
ent interplays arise between the two intervening processes depending on the initial 
densities B and Cq. As a result, unconventional and/or non-universal asymptotic 
behaviours may occur (B = C ), or the scaling regime may be absent (C > B ). 
In these models, the upper critical dimension is likely to be 4, as a consequence 
of the spatial structure induced by the catalytic reactants at late times: in 5-rich 
regions fragmentation (annihilation) will be frequent (absent) producing a high con- 
centration of A, whilst in C-rich regions the roles of these processes will be reversed 
producing voids of A. With simultaneous fragmentation and annihilation it seems 
difficult, however, to develop more quantitative arguments determining d c . In any 
case, new, unusual behaviours should also appear at lower dimensions in the pres- 
ence of catalysts. Even here, the mean-field description may prove adequate up to 
moderate times, before significant spatial fluctuations set in. 

Instead of the mean-field expressions (|33"D for the catalyst densities, one may 
use the corresponding ones in the diffusion-controlled regime {d < 4), which can be 
deduced, e.g., from f25| . In addition, other reactions, or time-dependences of the 
catalysts may be investigated. One may even replace one of the catalysts by an 
inhibitor. For these purposes one simply needs to use the results of Sec. 3.1. In 
practice, it may be that both linear and non-linear mass-loss processes operate. In 
such cases, models 1(C) and 2(C) would merge to yield a model with linear and 
non-linear annihilation terms. It should be possible to solve such a model using 
the methods of Sees. 2 and 3. It may also be that one has a source rather than a 
sink term, in which case the previous terms will have opposite signs leading to a 
competition between linear and non-linear contributions. 

Further directions of work are possible, in addition to those just mentioned. One 
may look for mass-dependent reaction-rates, the generic form for the annihilation 
terms then being: — Ja(x) dyf(x,y)a(y) for model 1, and —Ja(x)f(x) for model 
2. (In Sec. 2.1 we mentioned the case f(x,y) = (xy) a .) Another research avenue 
would be to investigate two species, say A and B, undergoing fragmentation and 
mutual annihilation. When A = B this model is identical to the single-species 
model of Sec. 2.1. With two species, a 'charge' conservation-law may also be present, 
e.g., A x + By — > A x _ y for x > y [B y - X for y > x) [§]. We have found |26| that in the 
absence of such a constraint both cluster densities A(t) and B(t) decay as 2/(Jt) 
(A = 1), independently of the initial masses A w and B w , i.e. there is no mass- 
difference residue. 

We expect the classes of models studied here, and their further developments, 
to have applications in various subjects, such as physics, chemistry, biology and, to 
some extent, population dynamics (see, e.g., [p3], |2"4]|). Model 1 (without catalysts) is 
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suitable to describe fragmenting systems with non-linear, reaction mass-loss. Models 
CI and C2, on the other hand, have the potential to describe a broad range of phe- 
nomena, with linear or non-linear mass-loss and with arbitrary time- dependences of 
the catalysts. We describe one particularly important generalization to these mod- 
els. In many processes of industry, heterogeneous chemistry, or biochemistry, such 
as food processing or food digestion (where the catalysts are enzymes), the efficiency 
of the reactions may depend on the available surface-area of the reacting fragments 



(see, e.g., P?fl ). It would be of great interest, therefore, to develop fragmentation 
models whose clusters are characterized not only by their volume (or mass), but 
also by their surface- area. From these, one could then build up fragmentation- 
annihilation models with surface-area dependent annihilation rates. An obvious 
difficulty in constructing such models is the dependence of the surface-area on the 
shape of the fragments. First steps towards this goal has been taken in very recent 



work on multivariable fragmentation [28 



It is interesting to note that there is some resemblance between fragmentation- 
annihilation processes and branching- annihilating random walks [p9fl . In the latter, 
particles diffuse and annihilate upon encounter and each particle gives birth to n 
offspring, at prescribed rates. A difference, though, is that the particles do not 
have a mass distribution. Depending on the branching and annihilating rates, these 
systems may evolve towards an empty state or an active steady-state. The existence 
of phase transitions in the long-time behaviour has been investigated for different 
values of n P^f . Depending on the densities of catalysts (which control the rates of 
the two processes involved), we have found that models CI and C2 can also exhibit 
quite opposite asymptotic behaviours: A — > for C > B , and A — > oo for B > C . 
Although we have not seen evidence for states characterized by steady values of the 
cluster- number density A(t), steady-state values were obtained for the cluster-mass 
density A\{t) (B > Co) and for the average cluster-mass L{t) (Co > B ). 

In conclusion, we hope that this paper will stimulate future work on fragmentation- 
annihilation processes, and believe that fruitful applications of the present, or related 
models (especially those with catalysts), can be found in various fields. 
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